Nematic and meta-nematic transitions in the iron pnictides 



S. Kasahara^'^, H.J. Shi^, K. Hashimoto^, S. Tonegawa^, Y. Mizukami^ , T. Shibauchi^, 
K. Sugimoto^'"^, T. Fukuda^'^''^, T. Terashima^, AndriyH. Nevidomskyy^ & Y. Matsuda^ 
^Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
'^Research Center for Low Temperature and Materials Sciences, Kyoto University, Kyoto 606-8501, Japan 
^Research & Utilization Division, JASRI SPring-8, Sayo, Hyogo 679-5198, Japan 
"^Structural Materials Science Laboratory, RIKEN SPring-8, Sayo, Hyogo 679-5148, Japan 
^Quantum Beam Science Directorate, JAEA SPring-8, Sayo, Hyogo 679-5148, Japan 
^Materials Dynamics Laboratory, RIKEN SPring-8, Sayo, Hyogo 679-5148, Japan 
^ JST, Transformative Research- Project on Iron Pnictides (TRIP), Chiyoda, Tokyo 102-0075, Japan 
^Department of Physics and Astronomy, Rice University, 6100 Main St., Houston, TX 77005, USA and 
* Present address: Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 

(Dated: July 5, 2012) 



Strongly interacting electrons can exhibit novel 
collective phases, among which the electronic 
nematic phases are perhaps the most surpris- 
ing as they spontaneously break rotational sym- 
metry of the underlying crystal lattice.^ The 
electron nematicity has been recently observed 
in the iron-pnictide^^ and cuprate*^^ high- 
temperature superconductors. Whether such a 
tendency of electrons to self-organise unidirec- 
tionally has a common feature in these super- 
conductors is, however, a highly controversial is- 
sue. In the cuprates, the nematicity has been 
suggested as a possible source of the pseudo- 
gap phase^^Ti^ whilst in the iron-pnictides, it has 
been commonly associated with the tetragonal- 
to-orthorhombic structural phase transition at 
Tg. Here, we provide the first thermodynamic 
evidence in BaFe2(Asi_a^Pa;)2 that the nematic- 
ity develops well above the structural transition 
and persists to the nonmagnetic superconducting 
regime, resulting in a new phase diagram strik- 
ingly similar to the pseudogap phase diagram 
in the cuprates^iiS Our highly sensitive mag- 
netic anisotropy measurements using microcan- 
tilever torque-magnetometry under in-plane field 
rotation reveal pronounced two-fold oscillations, 
which break the tetragonal symmetry. Combined 
with complementary high-resolution synchrotron 
X-ray and resistivity measurements, our results 
consistently identify two distinct temperatures — 
one at T*, signifying a true nematic transition, 
and the other at Ts{< T*), which we show to be 
not a true phase transition, but rather what we 
refer to as a "meta-nematic transition", in anal- 
ogy to the well-known metamagnetic transition 
in the theory of magnetism. Our observation of 
the extended nematic phase above the supercon- 
ducting dome establishes that the nematicity has 
primarily an electronic origin, inherent in the nor- 
mal state of high-temperature superconductors. 

In the iron pnictides, the antiferromagnetic transition 
is closely intertwined with the structural phase tran- 
sition from tetragonal (T) to orthorhombic (O) crys- 



tal symmetry. Although recent experiments, including 
neutron scattering^ ARPES^iii STM,^ and transport 
measurements^!^ have provided evidence for electronic 
anisotropy, these measurements were carried out either in 
the low- temperature orthorhombic phase ^i^iii where the 
crystal lattice structure has already broken C4 symme- 
try, or in the tetragonal phase under uniaxial strain^i^i^ 
that also breaks this symmetry. Therefore, the question 
remains open whether the electronic anisotropy can ex- 
ist above the structural transition without an external 
driving force, including under the superconducting (SC) 
dome. In the past, the nematic transition in the pnictides 
has been associated either with the orbital ordering^ — 
or with the spontaneous breaking of the Z2 Ising symme- 
try between two collinear magnetic ordering wave- vectors 
Q = (tt, 0) and (0, tt) J^ — Therefore determining the na- 
ture of the nematicity is a key to understanding the mi- 
croscopic origin of the lattice and magnetic transitions, as 
well as its possible connection with the high-temperature 
(high- Tc) superconductivity. ''^^ 

BaFe2As2-based materials are a prototypical family 
of iron-pnictides that play host to the superconductiv- 
ity upon ion substitution. Among them, the phase 
diagram of the isovalent pnictogen substituted system 
BaFe2(Asi_a,P^)2 (ref [H) is shown in Fig.[Tl The tran- 
sition to the antiferromagnetic ground state at Tn always 
coincides or is preceded by the T-0 structural transition 
at Tg. With increasing Tn decreases and goes to zero 
continuously at x = 0.30, indicating the presence of a 
quantum critical point (QCP).^^'^^ The superconduct- 
ing dome extends over a doping range 0.2 < x < 0.7, 
with maximum Tc = 31 K at the QCP*^ This system is 
very clean and homogeneous ^i^^i^S as demonstrated by 
the quantum oscillations observed over a wide x range 
even in the superconducting domeii^ This is important 
for the present work because impurities and sample in- 
homogeneities may otherwise wipe out the signatures of 
the phase transition. Thus BaFe2(Asi_a;Pa;)2 appears to 
be the most suitable system for studying the doping evo- 
lution of the intrinsic electronic and magnetic properties 
in the iron-pnictides. 

The magnetic torque r = jiioVM x ii" is a thermo- 
dynamic quantity, a differential of the free energy with 
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FIG. 1. Doping-temperature phase diagram of 

BaFe2(Asi_xPx)2. Solid circles, open triangles and open 
squares are the antiferromagnetic transition Tn , T-0 struc- 
tural transition Tg, and superconducting transition Tc tem- 
peratures determined by resistivity and thermal expansion 
measurements. Crosses indicate the nematic transition tem- 
perature T* determined by the present magnetic torque and 
synchrotron X-ray diffraction measurements. The insets il- 
lustrates the tetragonal FeAs/P layer. Xab = above T*, 
while Xab 7^ below T*, indicating the appearance of the 
nematicity along [110]t (Fe-Fe bond) direction. 



respect to angular displacement. Here V is the sample 
volume, M is the magnetization induced in the magnetic 
field and fio is the permeability of vacuum. Torque 
measurement detects the magnetic anisotropy with an 
extremely high sensitivity.^^ In particular, it provides a 
stringent test of nematicity. When H is rotated within 
the tetragonal ab plane (Fig. [2^, b), r is a periodic func- 
tion of double the azimuthal angle (f) measured from the 
tetragonal a axis (Fe-As/P direction, inset of Fig.[T]): 



^20 = ^fJ^oH'^Vlixo 



Xbb) sin 2^ - 2xab cos 2(j)] , (1) 



where the susceptibility tensor Xij is given by Mi = 
Xij Hj . In a system holding tetragonal symmetry, 
r2(f) should be zero because Xaa = Xbb and Xab = 0. 
Finite values of r2(f) appear if a new electronic or mag- 
netic state emerges that breaks the C4 tetragonal sym- 
metry. In such a case, rotational symmetry breaking is 
revealed by Xaa 7^ Xbb and/or Xab 7^ depending on the 
orthorhombicity direction. To avoid the formation of do- 
mains with different orientations in the ab-plane ('twin- 
ning'), we used very small single crystals (see Supplemen- 
tary Information). 

The upper panels of Fig.[2t depict the temperature 
evolution of the torque r(0) at fioH = 4T for opti- 
mally doped BaFe2(Aso.67Po.33)2 {Tc = 30 K). All torque 
curves are perfectly reversible with respect to the field 
rotation, r(^) can be decomposed as r{(j)) = + + 

H , where r2n(/) = ^2n0 sin2n(^ - (/)o) has 2n-fold 

symmetry with integer n. The middle and lower panels 
display the two- and four-fold components obtained from 



the Fourier analysis. The distinct twofold oscillations ap- 
pear at low temperatures, while they are absent at high 
temperatures. The fourfold oscillations (and higher- 
order terms) are observed at all temperatures (Fig.[2t, 
lower panels), but their amplitudes have negligible tem- 
perature dependence, indicating that they arise primarily 
from the nonlinear susceptibilities. 

As shown in Fig. [3^, the amplitude of the twofold os- 
cillation 1^201 is nearly zero at high temperatures and 
grows rapidly below T* :^ 85K in the optimally doped 
X = 0.33 crystals, followed by the saturation at lower 
temperatures. The anomaly at T* can also be seen by 
the synchrotron X-ray diffraction and in-plane resistivity 
measurements. Figure[3)3 depicts the temperature de- 
pendencies of the full width at half maximum (FWHM) 
and intensity of the Bragg peak below 200 K. Both the 
FWHM and the peak intensity change their slope in the 
vicinity of T*. Moreover, the temperature derivative of 
the resistivity dp/dT exhibits a peak at ~ T* (Fig.[3t). 
Here we stress that sharp synchrotron X-ray diffraction 
peaks, NMR spectra which show no indication of mag- 
netic ordering^^ and purely paramagnetic response above 
Tc (Fig. SI), all indicate that the sample inhomogeneity 
is highly unlikely to be an origin of the anomaly at T*. 

The above results clearly indicate that the tetragonal 
C4 symmetry, which is preserved at high temperatures, 
is broken below T*, demonstrating the formation of the 
electronic nematic phase at T*. The twofold oscillation 
below T* follows the functional form, r20 = ^2^ cos 20, 
meaning that Xaa = Xbb and Xab 7^ which indicates 
the nematicity along the tetragonal [110]-direction, i.e. 
Fe-Fe bond direction (Fig.[Tl inset). The X-ray FWHM 
at T < T* (Fig. [3)3, red circles) grows slightly larger than 
the linear extrapolation from above T* (dashed line), and 
is accompanied by the suppression of the X-ray peak in- 
tensity (green circles). This indicates a broadening of 
the Bragg peak below T*, implying that the nematicity 
to some extent couples to the orthorhombic lattice dis- 
tortion as discussed later. 

Figures [3]i-f show the temperature dependence of 
1^20 1, X-ray peaks and resistivity, respectively, of the 
underdoped material {x = 0.27). It undergoes the T-0 
transition at Tg ~ 60 K, and the magnetic transition 
well below Tc. |^20| is finite even at 200 K, and ini- 
tially increases with decreasing temperature, exhibiting 
a cusp-like peak at T* :^ 105 K, which we attribute to 
the nematic transition temperature for the following rea- 
sons. Analysis of more than five samples showed that 
the anomaly at 105 K (and also at Tg) is well repro- 
duced, while the torque curves above 105 K show some 
degree of variation depending on the sample (Fig. S2). In 
fact, when the temperature dependent part above 105 K 
is subtracted, 1^201 (Fig.[3li, open circles) exhibits very 
similar temperature dependence to that in x = 0.33 com- 
pound (Fig. [3^). In addition, it is apparent that the 
FWHM and intensity of X-ray peak (Fig. [3^) change their 
slope at T 105 K, quite similarly to those of the opti- 
mally doped compound at T* (Fig. [3)3). Moreover, the 
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FIG. 2. Torque magnetometry in the optimally doped x = 0.33 crystal, a, b, Schematic figures of the experimental 
configuration for the torque measurements under in-plane field rotation. The azimuthal angle (j) is defined as the field direction 
from the tetragonal a axis, c, Upper panels show the raw torque data r((/)) at several temperatures. Middle and lower panels 
are the two-fold r2cf) and four-fold components extracted from the Fourier analysis. 



hump structure can be seen in the dp/dT at around the 
same temperature (Fig.[3f). 

At T5 ~ 60K, 1^20! reaches a minimum and increases 
again at lower temperatures, indicating that 1^201 is sen- 
sitive to the T-0 structural transition. At temperatures 
above T*, the origin of the non-zero twofold signal is not 
clear. It may be due to the presence of local impurities or 
dislocations. In fact, local orbital ordering with C2 sym- 
metry around the impurity sites has been proposed. 

Figures[3^-i show the results for the parent compound 
X = 0, in which the structural transition and mag- 
netic order occurs at the same temperature T5 = T/v- 
Similar to underdoped compound, 1^20 1 is finite in the 
high temperature regime. As the temperature is low- 
ered, 1^20! increases gradually and then decreases with a 
cusp-like peak at ^ 170 K. At the same temperature, the 
synchrotron X-ray diffraction exhibits anomalies of the 
FWHM and peak intensity (Fig.[3]i). Thus, by the same 
reasoning as for the underdoped compound, the cusp-like 
behaviour of 1^20 |(T) is attributed to the nematic transi- 
tion; T* 170K. Behaviour of 1^20 |(T) at is also sim- 
ilar to underdoped compound but is more pronounced; 
\A2cf)\ is strikingly enhanced below T5, which is attributed 
to the larger distortion associated with the T-0 transi- 
tion. The resistivity shows a rather flat depenendece and 
no apparent anomaly is observed at T*. Both 1^201 and 
dp/dT exibit anomalies at around 50 K, which may in- 
dicate the presence of the magnetic or charge ordering 
deep inside the antiferromagnetic phase. 

We emphasize that in contrast to previous experi- 
ments, the present measurements were performed with- 
out applying external pressure or uniaxial stress. Thus 
our results provide thermodynamic evidence for the spon- 
taneous formation of the electronic nematic phase below 



T*, weh above the previously reported T-0 structural 
transition temperature Tg in this system. Moreover, the 
doping dependence of the nematic transition T*(x) dis- 
played in Fig. [1] indicates that nematicity develops inde- 
pendently of Tg {x) and persists over a wide range of dop- 
ing covering the non-magnetic superconducting regime. 

Clearly, there cannot be two nematic phase transitions 
at both Ts and T*, because the C4 rotational symme- 
try can only be broken once. The temperature T*(> Tg) 
marks the onset of the true phase transition, accompa- 
nied by the nematic two- fold torque component ^420 7^ 0. 
Then what happens to the structural transition at Tg? 
This question can be answered straightforwardly if one 
considers the Landau free energy expansion in terms of 
two parameters, the lattice distortion 5 and the nematic 
parameter ip ex A20, which can be written as follows: 

(2) 

with the terms in the first square bracket describing 
the first-order structural phase transition and the second 
bracket responsible for the (second-order) nematic phase 
transition. The temperature-dependent coefficients ts = 
{T-Ts^^^)/Ts^^^ Siudtp = {T-T^^^)/T^^^ were chosen such 
that in the absence of the coupling between the two or- 
der parameters, the structural transition occurs at lower 
temperature Ts^\< Tp^^). Thanks to the hnear couphng 
between the two order parameters, as expressed in the 
last term, both and S develop non-zero values below 
the transition temperature T* (Figs. [3^, b). On the other 
hand Ts ceases to be a true phase transition, since the 
C4 symmetry is broken on either side of T5, and lattice 
distortion S is non-zero over the entire temperature range 
(Fig. 13^). Instead, both 6 and i/j undergo a finite jump at 
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FIG. 3. Temperature dependence of the two- fold oscillation amplitude \A2(j)\ of the torque (blue circles), the 
FWHM (red circles) and the peak intensity (green circles) of the synchrotron X-ray Bragg reflection, and 
the in-plane resistivity p (black lines), a-c, Results for the crystal of x = 0.33, d-f, x = 0.27, g-i, x = 0. The nematic 
transition temperature T* and the met a- nematic transition temperature Ts are defined by arrows and vertical dashed lines. 
For underdoped (x = 0.27) and parent (x = 0) crystals, also shown are the corrected \A2cf)\ data (open circles in d, g) with 
subtraction of the smooth background of the two-fold oscillations above T* (Fig. S2). The X-ray are analyzed for [10,0,0]t (b) 
and [14,0,0] o (e, h) Bragg peaks and the x = data (h) shows a clear splitting below Ts (Fig. S3). The thin purple lines in c, 
f, i are the temperature derivative of resistivity dp/dT{T). 



Ts, as illustrated in Fig. SI We call this a meta-nematic 
transition, in analogy to the met a- magnetic transition in 
the theory of magnetism, where the magnetization un- 
dergoes a jump as a function of temperature or applied 
magnetic field, but remains non-zero on both sides of the 
transition. The analysis of the free energy shows that 
it exhibits a maximum at = and a single minimum 
at finite ip, as in the 2"^-order Landau phase transition 
(Fig. S4). Note that because of the coupling between 
the order parameters, both Ts and T* are renormalized 
compared to their initial values Ts^^ and Tp^^ (see Sup- 
plementary Information for details). 

To quantify the lattice distortion S experimentally, we 
have analysed the X-ray data by using two-peak fit- 
ting (Fig. S3), which reveals that the data in the region 
Tg < T < T* can be fitted with very small but finite 
S. The obtained results of S{T) can be reasonably re- 
produced within the framework of Eq. (|2]), and the same 
set of Landau parameters also fits well the temperature 
dependence of i/j (x (Fig.H]). We have thus estab- 
lished that the true thermodynamic transition occurs at 
T = T*, and is accompanied by the development of the 
non-zero values of both the nematic order parameter ip 
and the lattice distortion S. We note that similarly small 
but non-zero values of S have been recently reported in 
powder diffraction measurements of SmFeA(Oi_a;Fa;)j^ 



Note that this explanation is very generic and does not 
depend on the precise microscopic nature of the nematic 
transition, be it caused by Z2 spin-nematic ordering^ — 
or by orbital ordering. ^^^^ In the spin-nematic approach, 
the nematic instability is driven by thermal spin fluctua- 
tions above the SDW ordered phase, which would apply 
to the X < 0.30 regime in BaFe2(Asi_ajPaj)2- Such fluc- 
tuations have been detected by NMR3^ However, the 
fact that the nematic transition at T* occurs even for 
superconducting samples well above optimal doping, and 
far away from the SDW phase, provides a strong indica- 
tion that the nematic transition is unlikely to be asso- 
ciated with thermal fluctuations above long-range mag- 
netic order. On the other hand, in the case of orbital 
ordering, the nematic transition naturally occurs as a re- 
sult of polarisation between the Fe dxz and dyz orbitals, 
i/j oc {uxz — Tiyz), as supported by the recent ARPES^ 
and quadrupolar resonance measurements.^^ 

A large anisotropy in resistivity {pa — Pb)/{Pa + Pb) 
has been recently observed under the uniaxial stress in 
Ba(Fei_a^Coaj)2As2 at temperatures higher than Ts and 
even above the superconducting dome. In the light of 
the present results, this anisotropy can be explained as 
associated with the nematic transition at T*, rather than 
with a proximity to the structural phase transition. The 
application of a uniaxial stress will enhance the nematic- 
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ity and is likely to shift T* to a higher temperature. 



There is growing body of evidence that entanglement 
of the spin and orbital degrees of freedom leads to 
emergent novel electronic phases in the iron pnictides. 
The present temperature-doping phase diagram bears 
striking resemblance to that of high- Tc cuprates, in that 
the suppression of the antiferromagnetic ground state 
leads to the emergence of high- Tc superconductivity 
and electron nematic instability occurs well above the 
magnetic and superconducting transitions. Recent 
infrared studies of charge dynamics report the formation 
of a pseudogap in the excitation spectrum of optimally 
doped BaFe2(Asi_a;P^)2 below lOOK (Moon, S.J. et 
al, unpublished results). It is therefore likely that the 
nematic transition is related to the pseudogap formation, 
similar to the underdoped cuprates. These electronic 
properties may capture a universal feature essential for 
the occurrence of high- Tc superconductivity 
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Supplementary Information 

I. METHODS 

High-quality single crystals are grown by the self-flux 
method. To measure the magnetic anisotropy accu- 
rately, the microcantilever torque magnetometry is used 
(Fig. 2b). To avoid the torque oscillation arising from the 
misalignment, the magnetic field H is precisely applied 
in the ah plane within an error less than 0.1° by con- 
trolling two superconducting magnets and the rotating 
stage. 

In a nematic state, the domain formation with different 
preferred directions in the ah plane may occur as a conse- 
quence of the degeneracy in tetragonal crystal structure 
("twinning"). If the domain size is much smaller than 
the crystal size, the amplitude of T2(f) would be signifi- 
cantly diminished due to the cancellation of the two-fold 
oscillations with opposite sign arising from different do- 
mains. To avoid such an effect, we used very small single 
crystals with typical size ~ 70 x 70 x 30 /im^. 

In order to examine the detailed structural change of 
the underlying crystal lattice, high-resolution structure 
analysis by using the synchrotron X-ray crystallography 
techniques were performed on the same crystals used in 
the torque measurements and the high angle diffraction 
spots at tetragonal [10,0,0]t or orthorhombic [14,0,0]o 
were analysed. The sample temperature was controlled 
by an open flow cryocooler, whose temperature reading 
may have a difference from the actual sample tempera- 
ture by up to ^ 2 K, but this is sufficiently small for the 
purpose of discussion of the phase diagram in Fig.[Tl 

II. OUT-OF-PLANE MAGNETIC RESPONSE 

All of the single crystals used in the present study ex- 
hibit purely paramagnetic response above Tn and Tc, 
which was carefully checked by measuring the angular 
variation of the torque in H rotating within the ac plane. 
This excludes the possibility of extraneous magnetic im- 
purities in the samples. 

Figure Sla depicts the torque measured in field H ro- 
tated within the ac plane for underdoped {x = 0.27) 
single crystal at iiqH=4: T. The curves in this geome- 
try are perfectly sinusoidal, well fitted with r(T, 0) = 
A2e{T, H) sm26^ where A2e is the amplitude and 6 is 
the polar angle. The hysteresis component is less than 
0.01% of the total torque, indicating no detectable fer- 
romagnetic impurities. In this geometry, the difference 
Axca between the c axis and in-plane susceptibility yields 
a two- fold oscillation from T2e{0^T^H) with respect to 
rotation 

r2^ = i/ioi^VAxcaSm2^. (S-1) 

The iJ-linear dependence of A2e{H^T) / ji^H has null 
y intesect (Fig. Sib), which indicates a field-independent 



magnetic susceptibility, that is a purely paramagnetic 
response. This also reinforces the absence of ferromag- 
netic impurities. As the temperature is lowered, Axca = 
Xcc — Xaa increases gradually and then decreases with a 
peak at the T-0 structural transition Tg. 



HI. IN-PLANE MAGNETIC RESPONSE AND 
SAMPLE DEPENDENCE 

Figure S2 depicts the sample dependence of for 
underdoped {x = 0.27) compound. In both crystals, 
anomaly at T* is clearly seen. Solid symbols indicate 
the data, in which the temperature dependent ^420 at 
T > T* is subtracted by assuming linear-T dependence. 
After subtracting the high-temperature ^20, which may 
be due to the local impurities or dislocations, A2(f){T) 
shows nearly identical temperature dependence. Such a 
dependence is consistent with the raw data of x = 0.33, 
in which the high-temperature term is absent. 



IV. SYNCHROTRON X-RAY DIFFRACTION 
AND CRYSTAL STRUCTURE ANALYSIS 

Synchrotron X-ray diffraction were perfomed on 
BL02B1 at SPring-8 on the same crystals used in the 
torque measurements. The photon energy of the incident 
X-ray was tuned at 17.7 keV. A large cylindrical imag- 
ing plate and Rapid- auto program (Rigaku Corp.) was 
used to obtain diffraction data^^ The sample tempera- 
ture was controlled by an open flow cryocooler, whose 
temperature reading may have a difference from the ac- 
tual sample temperature by up to ~ 2 K, but this is suf- 
ficiently small for the purpose of discussion of the phase 
diagram in Fig. 1 in the main text. 

Figure S3a depicts the temperature dependence of 
diffraction angle 20 at orthorhombic (14,0,0)o for parent 
{x = 0) and underdoped {x = 0.27) crystals. The Bragg 
peak is analyzed by the two Gaussian curves (Fig. S3c). 
Below T*, the data can be fitted with two peaks, from 
which the orthorhobic distortion 6 = is estimated 
(Fig. S3b). 

The increase of the FWHM with decreasing T below 
200 K down to T* (Fig. 3b, e, h) is unusual because this 
temperature dependence is opposite to what is expected 
from the thermal vibration of atoms. This may indicate 
the presence of rather strong nematic fluctuations well 
above T*. We also note that the FWHM is supressed 
when the SC transition sets in. This implies the cou- 
pling between the nematic order and superconductivity. 
The peak of dp/dT at T* (Fig. 3c, f) suggests that the 
formation of the nematicity affects the electronic trans- 
port properties, possibly due to scattering off the nematic 
domain boundaries. 
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^(^^g) A)^(T) T(K) 

Figure SI. (a)The torque curve TaciO) in H rotating within the ac plane. {h)A2e / IJioH plotted as a function of iiqH . (c) 
Temperature dependence of Axca- 




peratures, which will become renormalized due to the 
electron-lattice coupling, as shown below. 

The saddle-point solution of the Landau free energy 
can be found from the following two equations: 



Figure S2. Temperature dependence of A24) for two samples 
of X = 0.27. The open symbols are the raw data. At high 
temperatures, the A2(/,(T) data are fitted with the T- linear 
dependence A^^(T) and from each torque curve t(0) the two- 
fold oscillation term with A^'^ (T) amplitude is subtracted to 
extract the reanalysed A2ct){T) (solid symbols). 



V. LANDAU FREE ENERGY ANALYSIS 

The Landau free energy is expressed in Eq. (2) in the 
main text in terms of the nematic order parameter ijj and 
the structural orthorhombic distortion 5 = (a — 6) / (a + 6) 
as follows: 

The quadratic coefficients are proportional to the reduced 
temperatures of the respective transitions: 



T-T, 

rpiO) 



(0) 



rpiO) 



(0) 



(S-2) 



where Ti°^ and T^^^ have a meaning; of the structural 
and nematic transition temperatures, respectively, in the 
absense of the coupling between the two order param- 
eters {g = 0). These are the "bare" transition tem- 



gS = 2tpilj ^ 4w ip^ 



(S-3a) 
(S-3b) 



The last equation can be thought of as expressing the 
structural distortion in terms of the nematic order pa- 
rameter ip: S = S{tlj). The free energy can then be ex- 
pressed in terms of ip alone: F{ip) = and is 
plotted in Fig. S4a for two different temperatures below 
(above) T5, shown in sohd (dashed) line. In both cases, 
the free energy has a maximum at = and a single 
minimum at finite ip, so that it behaves as in the 2"^- 
order Landau phase transition. This is a general result, 
provided the quartic coefficient w is larger than \u\ (in 
the opposite limit, the coupling between the 
two order parameters can render the nematic transition 
to be first order, but this does not alter the rest of our 
conclusions). 

We observe that the saddle-point solution, obtained 
by solving the equation dF/dip = 0, changes discontinu- 
ously when temperature goes through Tg, as illustrated 
in Fig. S4b. Consequently, structural distortion S also 
sustains a jump at T5, marking the onset of the meta- 
nematic transition, as explained in the main text. This 
is not a true phase transition because both order param- 
eters tp and S have finite values on either side of T5, as 
illustrated in Fig. 4a, b in the main text. 

The true phase transition occurs at a higher temper- 
ature T*, which can be obtained from the saddle-point 
equations by noting that sufficiently close to T*, the ip^ 
term can be neglected in Eq. (S-3b), resulting in 



for <T < T*. 



It then follows from Eq. (S-2a) that 



lltnr. 
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Figure S3. Temperature dependence of 20 (a) and 6 = If^l (b) for x = and 0.27. The color in (a) indicates the Bragg peak 
intensity, (c) Typical diffraction data (red circles) near the [14 0]o Bragg peak are fitted (blue lines) with two overlapping 
peaks (red curves) for T > T* (upper panels), T* > T > Tg (middle panels), and Ts > T (lower panels). 



This last equation has a solution for T < T* ^ and de- 
noting tp in terms of the reduced temperature tp = 

(T — Tp^^)/Tp^^ as in Eq. (S-2), we can express T* as 
follows: 

T* ^ T^^^) + . (S-6) 

Tp^^ has a meaning of the nematic transition tempera- 
ture in the absence of coupling to the lattice {g = 0). 
Because of this coupling, the transition is shifted to a 
higher temperature T* > Tp^\ 

Note that in the absence of coupling to the lattice 
{g = 0) the second-order electronic nematic phase transi- 
tion at Ti^^ would have been completely decoupled from 



the first-order structural transition at Tg < Tc , as illus- 
trated with the dashed lines in Fig. S5. However since 
both the nematic and orthorhombic order parameters 
break the C4 rotational symmetry of the crystal, it is 
natural to expect a non- vanishing couphng between the 
two. As a result, the orthorhombic distortion S develops 
a non-zero value below T*, rendering the transition at Tg 
a meta-nematic^ instead of a true first-order structural 
transition. This is corroborated by our X-ray diffraction 
studies, which find that the Bragg peak starts broaden- 
ing when temperature is lowered below T* , illustrated in 
Fig. 3b, e, h in the main text. 

The sharpness of the meta-nematic transition at Ts 
depends on the strength of the coupling g between the 
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Figure S4. Landau free energy analysis, Eq. (2) in the main text, with parameters u = 4:,v — 10, w — 8, g = 0.5 and ts{T),tp(T) 
expressed as in Eq. (S-2) via the unrenormahzed transition temperatures Ti^-* = 40 K, Tp^-* = 80 K. a, free energy as 
a function of the nematic order parameter ip for two temperatures, T = 60K and 80 K, on either side of the meta-nematic 
transition at Ts ~ 70 K. b, the saddle-point solution of Eqs. (S-3) for the nematic order parameter ip{T). 
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Figure S5. Landau free energy analysis with the same param- 
eters as in Fig. S4, except u = 1 and v = 20. Plotted with 
solid lines are the saddle-point solutions of Eqs. (S-3) for a, 
the lattice distortion 6{T) and b, the nematic order param- 
eter ip{T). The dashed lines denote the respective solutions 
in the absence of a coupling between the order parameters 
= in Eq. 2). 



order parameters relative to the typical jump in (5 at Tg, 
AS. To illustrate the point, we plot in Fig. S5 the saddle- 
point solutions of the Landau free energy using a different 
set of expansion parameters. As can be seen in Fig. S5, 
the meta-nematic transition becomes much smoother and 
less pronounced than in Fig. 4 in the main text, match- 
ing qualitatively the experimental data in the x = 0.27 
underdoped sample, see Fig. 4 in the main text. 

The above arguments are very general and do not de- 
pend on the microscopic origin of the electronic nematic 
order parameter i/j. It has been arguedi^riS that the ne- 
maticity may be driven by orbital physics, due to spon- 
taneous imbalance between the population of the Fe dxz 



and dyz orbitals, in which case oc (n^ 



This 



point of view is supported by recent angle-resolved pho- 
toemission^ and nuclear quadrupolar resonance (NQR) 
measurement^ and corroborated by ab initio calcula- 
tions Alternatively, the electron nematicity may be 
associated with the Z2-Ising spin orderingi^i^iS^ (some- 
times referred to as "spin-nematic"), in which case the 
order parameter can be written as oc {niA • ms) in 
terms of the sublattice magnetizations nij on the square 
lattice of Fe ions. In both cases, the order parameter ip 
is expected to couple linearly to the orthorhombic lattice 
distortion S. 

While our free energy analysis holds independently of 
the origin of nematicity, the fact that the two-fold os- 
cillations of the torque magnetization are observable in 
the strongly overdoped region {x = 0.5), far away from 
antiferromagnetic region in the phase diagram (see Fig. 1 
in the main text) suggests that this nematic signal is un- 
likely of magnetic origin and that orbital physics likely 
plays an important role. 
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